Dephasing of Si spin qubits due to charge noise 
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Spin qubits in Silicon quantum dots can have long coherence times, yet their manipulation relies 
on the exchange interaction, through which charge noise can induce decoherence. Charge traps near 
the interface of a Si heterostructure lead to fluctuations in the quantum-dot confinement and barrier 
potentials, which cause gating errors and two-spin dephasing. We quantify these effects in Si double 
quantum dots using a realistic model of noise. Specifically, we consider both random telegraph 
noise from a few traps (good for dots grown on submicron wafers) and 1// noise from many traps 
(good for larger wafers appropriate for quantum dot arrays). We give estimates of gate errors for 
single-spin qubit architectures and dephasing in singlet-triplet qubits. 



Solid-state spin-based qubits [T] are believed to be 
promising for scalable quantum computation (QC) . Spins 
are weakly coupled to the environment leading to long 
spin coherence times, while semiconductor-based struc- 
tures, such as quantum dots (QDs), can in principle be 
straightforwardly scaled up. Among semiconducting host 
materials. Si stands out with small spin-orbit interac- 
tion, no piezoelectric electron-phonon coupling, and weak 
(which could be further reduced through isotopic purifi- 
cation) hyperfine coupling, so that single electron spins 
have excellent quantum coherence properties, with a co- 
herence time for donor electron spins in Si:P of up to 
300 ms [5]. This coherence time is the longest among 
sohd state qubits, and considerably exceeds the electron 
spin coherence time of ~ 1/US in GaAs QDs [21 14], which 
is plagued by hyperfine interaction with nuclear spins. 
Moreover, the sophisticated Si fabrication technology is 
a potential asset in scaling up Si-based QC architectures. 

At present, lateral quantum dots have been fabricated 
in Si structures such as Si/SiGe quantum wells [5] and 
Si/Si02 heterojunctions [B]. An important experimental 
and theoretical issue concerns the spin coherence prop- 
erties in a single or double QD near a Si/Si02 interface, 
where the presence of defects such as Pb centers is in- 
evitable and the oxide is well known as a source of charge 
noise [7]. Will the superior spin coherence properties in 
bulk Si hold up? Even if charge fluctuations do not af- 
fect single spins significantly, when spins are coupled via 
exchange interaction, they can be adversely affected by 
charge noise because exchange is Coulombic in nature [H] . 
Since the main perceived advantage of Si spin qubits is 
the expected long spin coherence time because of weak 
hyperfine coupling (and the potential for isotopic purifi- 
cation by removing Si-29 nuclei, thus further reducing 
the coupling to nuclear spins), it is crucial to consider 
the decoherence of Si spin qubits due to the electrostatic 
noise or fluctuations in the environment affecting the ex- 
change coupling in the system, which is independent of 
nuclear spin considerations. 

In this Letter we quantify the effect of charge noise 
on exchange-coupled spin qubits in Si double quantum 



dots (DQDs). In particular we calculate two-qubit gate 
errors when single spins are used as qubits [Tj, and de- 
coherence rates when two-spin singlet and triplet states 
serve as encoded qubit states [3]. Below we first give a 
brief discussion of charge noise near a Si/Si02 interface, 
then set up the DQD Hamiltonian. We calculate the 
exchange coupling in the Heitler-London approximation 
and determine the modifications to the exchange due to 
variations in the barrier and level detuning. Finally we 
discuss dephasing rates in singlet-triplet qubits. 

In a Si MOSFET structure dangling bonds {Pb cen- 
ters) near the interface act as charge traps, capturing 
electrons from a nearby source such as a 2-dimensional 
electron gas (2DEG) and re-emitting them. Charge traps 
reside in the oxide layer up to a few nm from the inter- 
face, within tunneling distance from the 2DEG. A trap 
is characterized by an activation energy Et, which is the 
energy required to change the charging state of the trap. 
For example, Pb centers generally lie below the Fermi en- 
ergy [5J [TU] , thus the traps are mostly occupied, so that 
Et — Ep, where Ep is the Fermi energy of the 2DEG, 
corresponds to the energy the trap needs to emit an 
electron to the Fermi surface of the 2DEG. The mean 
capture and emission times (tc) and (rg) of the traps 
satisfy {tc) / {re) = e^^'-^'')/'^^^. These times depend 
sensitively on temperature and gate voltage, and in a Si 
MOSFET at low temperatures they can range from <1 
ms to >1 s The number of carriers in the trap typ- 
ically fluctuates by one. Charge fluctuations due to one 
trap are represented by a random function l^V(t). We 
take Ay to have the same magnitude for all traps and ne- 
glect interactions between traps. In a classical transistor, 
the alternate capture and emission of carriers by individ- 
ual defect sites generates discrete switching events in the 
current through the device called random telegraph noise 
(RTN). The spectral density of RTN is 
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In submicron devices there are typically only a few 
traps [12] • Yet scalable quantum computation requires 
the fabrication of large QD arrays. In such devices, which 
can be ^ l/im in length scale, the ensemble of traps has a 
distribution of activation energies Et, and thus have dif- 
ferent capture and emission times. Such a sum of many 
random telegraph signals yields a 1// spectral density 
[H [H], a fact well documented in MOSFETs [13 IH] 
and other microelectronic devices [THJ \T7\. For exam- 
ple in Ref. [18], where 1// noise was observed, the trap 
density ~100 - 1000/im^^, so that across a 100 nm sec- 
tion one expects a few fluctuators. This is the case in 
small samples where discrete switchings are resolved [T^ . 
For 1// noise, as its name indicates, the spectral density 
takes the form 5*^(0;) = aksT/uj, where a represents the 
strength of this noise and a cutoff loq is introduced at the 
small-o-i limit, often taken as the inverse of the measur- 
ing time. In the present study we consider both RTN 
and 1// noises. 

To explore how charge noise affects exchange-coupled 
spin qubits in Si, we study a DQD in a 2DEG in a 
MOSFET structure grown along the z direction. The 
right and left dots are located at Rr^l = (i-'^OjO, 0) 
respectively. The two-electron Hamiltonian is Hq = 

( T,^=l,2 ^^'^ + ^Q ^) + ^ee, whcrc T IS the kinetic energy 
and Vq the confinement potential. 
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with a and b the Fock-Darwin radii for in- and out-of- 
plane confinement, and and the in- and out-of- 
plane Si effective masses. Vx is given below for symmetric 
and asymmetric dots. The Coulomb potential between 
electrons at ri and r2 is T4e — e^/(e|ri — r2|), where 
e = (esi + esiOa )/2, accounting for the effect of the image 
charge in the oxide. 

In a DQD charge fluctuations can cause the height 
of the potential barrier to fluctuate and introduce ran- 
dom offsets between the bottoms of the two dots. For 
single-spin qubits, the exchange coupling J controls in- 
teractions between qubits and is pulsed on and off by 
a gate. Here we choose a confinement potential = 

2m7c? ( ^ 2X^° ) ^ ' '^hich gives rise to a barrier potential 
Vbi — {h^ / Smttt'^) . In the Heitler-London approxima- 
tion 
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where I{R) = /„ dqe " T ^, 
with Jo the zeroth-order Bessel function and Erfc the 
complementary error function. 

In general the total electrostatic potential that defines 
a quantum dot contains contributions from all possi- 
ble sources - gates, defects, static dopants, etc. Each 



source contributes to the barrier as well. The exchange 
J is a function of all the contributions and can be writ- 
ten as J{Vbi,Vb2, ■■■)■ The total fluctuation in J is ap- 
proximately A J — J2j{9J/dVbj) AVbj, where j labels 
the different barrier contributions. The qubit's sensitiv- 
ity to fluctuations in barrier j is therefore encapsulated 
in dJ/dVbj. Here we consider only two contributions 
to the barrier - from the original confinement potential 
Vbi = Vx and from the trap potential Vb2 = Vr- While 
Vx depends on device parameters, Vt is determined by 
the random distribution of traps. When the trap fluctu- 
ates only Vb2 changes, thus AJ « {dJ/dVb2) AVf,2- When 
the potential of the gate electrodes fluctuates to give a 
nonzero AVf,i, it causes a fluctuation in the exchange of 
A J « {dJ/dVbi) AVbi. With Vbi and Vb2 having different 
geometry and dynamics, and exchange J depends nonlin- 
early on these variables, generally dJ/dVb2 ^ dJ/dVbi- 

We now estimate the fluctuation A J for a double dot 
under the influence of RTN, since a fluctuating J leads 
to a gating error J J dt that results in unwanted entan- 
glement in spin states W- For a single DQD we expect 
1-4 traps [12], of which perhaps only one will be close 
to either dot. We consider a trap at Xt just outside to 
the left of the left dot, at a depth d in the substrate. 
In the Thomas-Fermi approximation, the screened trap 
potential is 
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where qxF = 2/0^ [19 , with ^ 3nm the effective 
Bohr radius in Si. In the Heitler-London approxima- 
tion we find J ^ J + AJ, where A J ~ (Au^ -t- Aw/j — 
2Aw)/sinh(2X^/a2), with /^vl,r = {L,R\Vt\L,R) and 
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{L\R) ' ^'^^ have used |L, R) for the left and 
right dot wave functions. We consider a QD with an 
electron state radius of 8nm. It is expected that a much 
larger area of radius ~50nm around the dot would be 
depleted so that any trap in this area will not be charge- 
active. For a realistic 50 nm interdot separation this 
gives us Xt — —75 nm. With d taken to be negligible 
and e w 8 we find J w 0.4/LteV and AJ/J « 7 x 10"^ 

Energy-level fiuctuations have been observed in GaAs 
QDs tl6j, in Si MOSFETs [El Hg [HI [201 [21] , and in Si 
nanowires [TT]. For a GaAs QD the level fluctuation is 
found to be ~0.07-0.16/j,eV [TB], while for a reasonably 
quiet Si MOSFET QD it is ~0.45AteV [2D]. However, 
the length scale of these fluctuations has not been de- 
termined at present. Simultaneous fluctuations in the 
dot potentials of the same magnitude do not affect the 
exchange. In other words, only short wave length fluc- 
tuations, which affect the interdot bias or barrier height, 
affect the value of J and cause dephasing. Neverthe- 
less, based on the similar values of level fluctuations in 
GaAs and Si QDs, we conclude that, while Si can have 
much better one-qubit coherence properties than GaAs, 
for two-qubit gates it will have similar gate errors due to 
charge fluctuations. 
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Figure 1: Dephasing of a singlet-triplet qubit in a Si/Si02 
DQD with radius a = 8.2nm and T =0.1K for RTN and 
1/f noise. The vertical axis shows the surviving amplitude 
Pi2 (t) I pi2 (0) as a function of time. For RTN we have chosen 
T = O.lms. For 1/f noise the surviving amplitude is shown 
for different values of dJ/dV, with Sv{lo) « 3 x W'^^Y'^ /lu 
estimated from Ref. [20]. The inset shows the exchange cou- 
pling J as a function of bias detuning for a DQD with a = 
8.2nm and interdot separation 50nm. 



A logical qubit can be encoded in the two-electron sin- 
glet and unpolarized triplet spin states in a biased regime 
for a DQD [3]. In this case, if a finite J is maintained to 
avoid singlet-triplet mixing by an inhomogeneous mag- 
netic field, charge fluctuations could lead to pure de- 
phasing between the qubit states. To assess this dephas- 
ing in Si, we calculate J with a tilted model potential 
14 = (hujo/^a^) Mm[{x - Xo)^, {x + Xq)^] - eEx, where 
the external electric field E raises the energy of the left 
dot with respect to the right dot. Here J is controlled by 
applying a bias potential energy e ~ 2eEXQ between the 
two dots and can be very large, as shown in the inset of 
Fig. [l] It is given approximately by J ~ 2P/e, where i 
is the tunnel coupling and e is the energy splitting of the 
two lowest singlet states. 

In a tilted potential, variations in the exchange are 
given by A J = (dJ/di) At + (dJ/de) Ae. The tunnel 
coupling is determined mostly by the interdot distance 
and details of the interdot barrier. It varies only slightly 
within the bias regime that we work with, and is treated 
as an independent variable from e. Therefore we can 
write A J = dJ/deAe. Ae has the same electrical origin 
as the level fluctuations AV and has contributions from 
both the electrodes and the charge traps: Ae = AVhi -I- 
AVb2- This linear relationship allows us to evaluate the 
effect of the charge traps based on the gate potentials 
since dJ/dVb2 = dJ/dVbi — dJ/de. If we focus only on 
noise from charge traps, A J — {dJ/de)AV. 

The singlet and triplet states define a 2D Hilbert space, 
in which the effective Hamiltonian can be written simply 
as Heff = [J-l- AJ{t)] az- In other words, charge fluctu- 
ations lead to pure dephasing in the singlet-triplet qubit 



and a decay in the off-diagonal element pi2 of the system 
density matrix. If the double dot qubit is only affected 
strongly by a few nearby charge traps, the corresponding 
charge noise can be treated as RTN, which is strongly 
non-Gaussian. We find that |22| 
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where rj = ^/{J + AJ)^/h? + {1/t^). Since (J -I- 
AJ)/h > l/r typically we have 77 sa ( J -f AJ)/h. Two 
time scales enter Eq. ([s]) - one associated with precession 
and one with decay. While the precession time scale is 
set by J and is fast (shorter than ns), the decay time 
scale is set by the switching time of the fluctuator. This 
switching time can be very long and can exceed the spin 
lifetime in Si. For example, experiments on a Si/Si02 
nanowire at dilution refrigerator temperatures [11] have 
found T ranging from 0.1 ms to hundreds of ms depending 
on the gate voltage. From this perspective working with 
small samples having few traps is obviously beneficial. 

A scalable quantum computer will no doubt be made 
on large wafers, which are usually subject to 1/f noise, 
typically Gaussian in MOSFETs . In this case we can 
write /Oi2(i)/pi2(0) = e~^^*\ where 
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For 1/f noise the low frequency part of the 
spectrum dominates dephasing, so that x(t) ~ 
{dJ/dVf{akBT/2h^)t^\a.uiot. The time evolution of 
the off-diagonal element of the density matrix is given 
by Fig. [1] where the electron temperature has been set 
to 100 mK, the typical value in a dilution refrigerator. 
Figure [T] shows that dephasing kicks in quite rapidly for 
the parameter regime we choose, ranging from 0.1 /iS to 
10 /IS. The dephasing time can be enhanced in three 
ways. One can go to ever lower temperatures to reduce 
Sv{oj), yet this is limited by available technology. One 
can also separate the dots further, since larger interdot 
distances lower dJ/de. The best approach may be to 
apply a magnetic field and find the so-called sweet spot, 
at which dJ/de — and the double dot is insensitive to 
charge noise |53]. Finding such sweet spots may, how- 
ever, be problematic since quantitatively exact numeri- 
cal estimates of J(V), taking into account all electrostatic 
contributions, would be difficult, if not impossible, in a 
scaled up architecture. 

In summary, we have studied fluctuations in the ex- 
change coupling and dephasing of spin qubits in Si DQDs. 
Single-spin qubit systems in Si are expected to have sim- 
ilar gating errors to GaAs. For singlet-triplet qubit sys- 
tems there is a significant difference in the two-spin de- 
phasing between submicron devices, in which RTN is ex- 
pected to be dominant, and large devices, in which 1/f 
noise dominates. 
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